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Abstract 

We study the long time behaviour of local observables following a quantum quench in 
1+1 dimensional conformal field theories possessing additional conserved charges besides 
the energy. We show that the expectation value of an arbitrary string of local observables 
supported on a finite interval exponentially approaches an equilibrium value. The equilibrium 
is characterized by a temperature and chemical potentials defined in terms of the quenched 
state. For an infinite number of commuting conserved charges, the equilibrium ensemble is 
a generalized Gibbs ensemble (GGE). We compute the thermalization rate in a systematic 
perturbation in the chemical potentials, using a new technique to sum over an infinite number 
of Feynman diagrams. The above technique also allows us to compute relaxation times for 
thermal Green’s functions in the presence of an arbitrary number of chemical potentials. In 
the context of a higher spin (hs[A]) holography, the partition function of the final equilibrium 
GGE is known to agree with that of a higher spin black hole. The thermalization rate from 
the CFT computed in our paper agrees with the quasinormal frequency of a scalar field in 
this black hole. 
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1 Introduction and Summary 

The study of thermalization in closed interacting quantum systems has a long history (see, 
e-g- [I] for a review). It has been known ever since the celebrated work of Fermi, Pasta and 
Ularn (FPU) that interacting classical systems need not necessarily equilibrate. The question 
of finding sufficient conditions for thermalization in quantum systems is also an open one. 
Recently, the advent of holography has linked the issue of thermalization in strongly coupled 
quantum held theories to another important, classical, problem of black hole formation (see, 
e.g. M and references therein). In the latter setting too, the issue of gravitational collapse 
of a given matter distribution is rather nontrivial; indeed there is an interesting debate in 
the current literature (see, e.g., 0) regarding the fate of perturbations in anti-de-Sitter 
spacetimes. 

In this paper, we will focus on two-dimensional conformal held theories (CFTs) on an 
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infinite line a E M. We will consider the system at t = 0 to be in a “quenched state’ Q 

OO 

iVb) = exp [~€ 2 H - j ~^e n W n \\Bd) (1) 

n =3 

Here | Bd) is a conformal boundary state; the exponential factors cut off the UV modes to 
make the state normalizable. W n denote the additional conserved charges in the theory^ 
This choice of the quenched state is a generalization of that in [5] for which e n = 0, for n > 3. 
The wavefunction for t > 0 is given by 

IV’W) = exp[-iHt]\ij)o) (2) 

The questions we will explore, and answer, are: what is the long time behaviour of 
various observables in \i/j{t))7 In particular, does the expectation value of an operator (or 
a string of operators) approach a constant? If so, (i) is the constant value characterized by 
a thermodynamic equilibrium, and (ii) what is the rate of approach to the constant value? 
More generally, we would also address, and partially answer, the questions: how does the 
existence and rate of thermalization depend on the initial state and the choice of observables? 


Thermalization We find in this paper that the expectation values of local observables 
(supported on a finite interval A : a E [—1/2,1/2]) asymptotically approach (see (fT2jl for the 
precise statement) their values in an equilibrium ensemble, 

Peqm = 77 exp [~/3H - ^ fi n W n ], Z = Tr exp [~/3H - Y" n n W n ] ( 3 ) 

n n 

whose temperature and chemical potentials are related to the cutoff scales in (J5]) as follows 

(3 = 4 : 62 , p n = 4e n , n = 3 , 4 ,... ( 4 ) 

The relations ((4]) are uniquely dictated by the requirement that the expectation values of 
the conserved charges H, W 3 , W 4 , ... in the initial state match those in the mixed state ()3j) 
(see (l52]l ). In the absence of the extra parameters e n , n = 3 , 4 , ... this result is derived by the 
elegant method of conformal transformations [5]. In the presence of these parameters, this 
method is not available; in this paper, we deal with the extra exponential factors in terms 
of an infinite series and do a resummation. 

4 In the original sense of the term, a quantum quench is defined as a sudden change from a hamiltonian 
Hq to a hamiltonian H which governs further evolution for t > 0. The system is assumed to be in the ground 
state of Hq at t = 0, which serves as an initial state for subsequent dynamics; the dynamics is nontrivial 
since the initial state prepared this way is not an eigenstate of H. In this paper, as in [2, we will mean 
by a “quenched state” simply a pure state which is not an eigenstate of the Hamiltonian H . The kind of 
quenched state defined in 0 is sometimes said to describe a global quench or a homogeneous quench, as the 
state is translationally invariant. We will briefly mention inhomogeneous and local quenches in Section [G] 
5 For the purposes of this paper, we will identify them with lT n -charges of 2D CFT, n = 3,4,... (with 
W 2 = H), although much of what we say will go through independent of this specific choice as long as these 
charges mutually commute and are defined from currents which are quasiprimary fields of the conformal 
algebra. 
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We emphasize that the thermalization we found above persists even when we have an 
integrable model with an infinite number of conserved charges. Relaxation in integrable 
systems has been found in recent years in the context of, e.g., (a) one-dimensional hardcore 
bosons | 6 j, (b) transverse held Ising model [7], and (c) matrix quantum mechanics models 
| 8 j. The equilibrium ensembles in this context have been called a generalized Gibbs ensemble 
(GGE). Our present result on integrable conformal held theories adds to the list of these ex¬ 
amples. Interestingly, the thermalization we find works even for free conformal held theories, 
e.g. a free scalar held theory]^ 

With the above identihcation of parameters, we will rewrite the initial quenched state 
© henceforth as 


|Vh) = exp [-{PH - ^2fi n W n )/4:]\Bd) (5) 

n=3 

We find the following specihc results: 

1. Thermalization time scale for single local observables: We find that at large 
times 


= Tr (0 fc (O)p egm (^,/ii)) +a k e lkt + ... ( 6 ) 

where </>fc(cr) is an arbitrary quasiprimary held (labelled by an index k ). Below we compute 
the thermalization exponent 7 ^. in a perturbation in the chemical potentials and to linear 
order it is given by 



A k + p nQn,k + 0 (jj 2 ) 

n 


pn 


Pn 


(7) 


Here A k = h k + h k is the scaling dimension and Q ntk are the (shifted) BG-charges (see (IHTD 
for the full definition) of the held <f k (in case of primary helds) or of the minimum-dimension 
held which appears in the conformal transformation of <f k . To obtain this result, we perform 
the infinite resummation mentioned below (J4J). At large times, the perturbation series for the 
one-point function in the chemical potentials exponentiates (see (137| ). to give the corrected 
exponent in the above equation. In various related contexts, finite orders of perturbation 
terms in chemical potentials have been computed before ie dam. Our finding in this paper 
is that at large times, there is a regularity among the various orders leading to an exponential 
function as in (see Section Y2.2.2\ for details). 

Universality: In the case of zero chemical potentials, it has been noted in ra. that although 
the relaxation time 77 = 7 ^ 2 / ( 2 A k ) = 2n/3/(A k ) is non-universal (in the sense that it depends 
on the specihc initial state (JTj) ), the ratio of relaxation times for two different helds, namely, 
rfci/rfc 2 = A k 2 /A kl is universal (it depends only on the CFT data and not on the initial state 
and is hence expected to be valid for a general class of initial states). In the presence of the 

6 This happens essentially due to the fact that we consider here thermalization of local observables and 
that local held modes are mutually coupled even in a free held theory. Thermalization happens at times 
greater than the scale of localization, as we will see below. 
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additional cut-off parameters e^i — 3,... in the initial state ([I]), the ratio r^/r^ = r jk 2 / r 7k 1 = 
(A fc2 + Y.nl 1 nQnte)/(^k 1 + J2nPnQnM) is > however, not independent of the initial state. 

However, as we will briefly discuss in Section [ 6 j for a large class of quench states 
(e.g. where the energy density is uniform outside of a domain of compact support) the 
/3-dependence of 17 , in the absence of chemical potentials, can be understood as the depen¬ 
dence on the uniform energy density (see a related discussion in The time scales 77 , 

therefore, do have a limited form of universality in the sense that it depends on a rather 
robust feature of the initial state. Our calculations in this paper leads us to believe that this 
feature will continue in the presence of chemical potentials, in the sense that the additional 
dependence of the time scales 1 / 7 *, on the p n is fixed by the charge densities corresponding 
to the additional conserved charges. We hope to address this in p4 j. 


2. Multiple local observables, reduced density matrix: Besides the one-point func¬ 
tions discussed above, it turns our that we can demonstrate thermalization of all operators 
in an interval A of length l. It is convenient to define a ‘thermalization function’ /^(t) [32j 
as 


T (+\ — tW " /yr (R 'O - Tr (pdynA^P eqm,A (/3, pn)) 

lA\t) ~ i-APdyn,A\ t )Peqm,AyP^n)) ~ “~ 77\2YV~( / a \2X\V 2 

[TltpdynAt) )*- V {Peqm,A{Pi Pi) )\ 

Pdyn,A(f) Wa |"0(t)) ^"0(t) |, Peqm,A (/3, Pn) Wj4 PeqmifliPi) 


( 8 ) 


Here p = p/ \/Tr p 2 denotes a ‘square-normalized’ density matrix FFI We show below that at 
large times the thermalization function has the form 


I A (t) = 1 - a(!) e~ 2lmt + ..., I = 1/(3 (9) 

where refers to the exponent ((7| for the operator (f) m with minimum scaling dimension^ 
a{l) is computed as a power series in l which we find using the short interval expansion, 
valid for l • Cl, i.e. I <C (3. 

Two immediate consequences of (J9J) are 

(i) Thermalization of an arbitrary string of operators: Note, from (J9J) , that 


iAit) 


t—> 00 _ 

- > 1 , 


( 10 ) 


Since the square-normalized density matrices can be regarded as unit vectors (in the 
sense of footnote ED, and I Ait) can be regarded as the scalar product Pdyn,A{t)'Peqm,A, 
(flOli clearly implies 


/s / 4 -\ ^^ 

Pdyn,A\t) ^ Peqm,A 


(ii) 


7 Note that operators in a Hilbert space H can themselves be regarded as vectors in H x H*; under this 
interpretation Tr(T B) defines a positive definite scalar product. With this understanding, we will regard 
the hatted density matrices as unit vectors. 

8 Throughout this paper, we will consider field theories with an infinite spatial extent. The entire Hilbert 
space is assumed to be of the form H 4 ® H 4 . Tr^ implies tracing over H^. 

9 We will assume here that the spectrum of such A’s is bounded below by a finite positive number. In 
case of a free scalar field theory, we can achieve this by considering a compactified target space. 
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This implies the following statement of thermalization for an arbitrary string of local 
operators (with a±,a 2 ,... G A) 


(^(t)|0(oi, t 1 )0(a 2 , t 2 ). .#(<)) = Tr(p dvn)A (t)0(<r 1 ,t 1 )0(a 2 , h)-~) 

Tr(p eqmtA 0 (ai,ti) 0 (a 2 ,t 2 )...). (12) 


(ii) Long time behaviour of reduced density matrix: 

Carrying on with the interpretation of 7^(t) as a scalar product, we can infer following 
asymptotic behaviour of Pdyn(t) from (JSJ): 

Pdyn,A(t) = Peqm,A((3i Pi) — O' e 2lrnt + ...)+ Q ^\/2a e 7m4 + ...^ (13) 

where Tr(Q 2 ) = 1, T^(Peqm,A(/3, Pi)Q) = 0. We will specify further properties of Q 
later on. 

Importance of local observables: In case of a free massless scalar field, it is easy 
to show that quantities like (^>{t)\a\o\\x/{t)) perpetually oscillate and never reach a 
constant (see a related calculation in [8]). The modes a n represent Fourier modes and 
are non-local. Indeed, as mmm showed, in the absence of chemical potentials, 
the exponential term in ([2]) is e ~ 2 'irn(t-i/ 2 ) an q iq ie thermalization sets in only after t 
exceeds 1/2. Thus, for l = oo, there is no thermalization, which is consistent with 
the above observation about perpetual oscillations. We expect the form e ^ 2 7m(f-F 2 ) to 
continue to hold in the presence of chemical potential^, since the effect of the chemical 
potentials on the exponent y*, can be viewed as a shift of the anomalous dimension 
5A k = PnQn,k + 0 (p 2 ) (see, e.g. (JMD). This shows that, as in the case of zero 
chemical potentials, equilibration sets in only after t exceeds 1/2. We will see a similar 
phenomena next in the context of a decay of perturbations to a thermal state. 

3. Decay of perturbations to a thermal state: We compute (see Section H] for details) 
the time-dependent two-point Green’s function G + (t,l; /3, p) for two points spatially sepa¬ 
rated by a distance l. We find that for — l /3, the time-dependence is exponential, 
with the same exponent as in (J6]l: 

G+(t, l, /3, p) = -^Tr (0 fc (Z,t)0fc(0,0)e" /3i/_E " MnW "’ 1 ) = const e~ lkt (14) 

Note that the above thermalization sets in for t > l. For t < l, the two-point function has 
an exponential decay in the spatial separation (see Section [4] and Figure |3|). 

The computation of the above relaxation times in the presence of an arbitrary number of 
chemical potentials uses the technique, described above, of summing over an infinite number 
of Feynman diagrams, and is one of the main results of our paper. 

10 Although, in the short-interval expansion employed in this paper to derive (fUjl. which uses t gg/3 l, 
such an ^-dependence in the exponent cannot be easily seen from the pre-factor a(l) unless one sums over 
an infinite orders in l. 


5 





4. Collapse to higher spin black holes: In [T7] the bulk dual to the time-dependent 
state ([2]) corresponding to initial condition ([5]), for large central charges, has been constructed 
in the case of zero chemical potentials. The dual geometry corresponds to one half of the 
eternal BTZ (black string) geometry, whose boundary represents an end-of-the-world brane. 
In [18] the result has been extended to the case of non-zero angular momentum and a 
Chern-Simons charge. In case of an infinite number of chemical potentials, a bulk dual to 
the equilibrium ensemble (E]) has been identified, in the context of the Gaberdiel-Gopakumar 
hs(A) theory [H5], as a higher spin black hole with those chemical potentials |2Uj . It is natural 
to conjecture mm that the time-development (J2]) should be dual to a collapse to this higher 
spin black hole. At late times, therefore, the thernralization exponent found above should 
correspond to the quasinormal frequency of the higher spin black hole. We find that (see 
Section 0 and TO this is indeed borne out in a specific example. 

The plan of the paper is as follows. The results 1, 2, 3 and 4 above are described in 
Sections El H El and El respectively. The resummation of an infinite number of Feynman 
diagrams (corresponding to insertions of arbitrary number of chemical potential terms) is 
discussed in Section 12 .2.21 which uses results in Appendix El The calculation of the overlap 
of reduced density matrices in Section E] needs the use of the short-interval expansion, which 
is described in Appendix [B] In Section [6] we present our conclusions and make some remarks 
on inhomogeneous quench m 

2 One-point functions 

In this section we will consider the behaviour of the following one-point functions of a 
quasiprimary held (f) k {a) 


(4* h eqm = Tl (0fc(o")p e g( m (/3, /X n )) (15) 

We will briefly recall how these are computed in the absence of the chemical potentials [5j ‘22] . 
The first expectation value corresponds to the one-point function on a strip geometry, with 
complex coordinate w = a + ir, a 6 (— 00 , 00 ), r G (—/l/4,/3/4) where r is eventually to 
be analytically continued to r = it. This can be conformally transformed to an upper half 
plane by using the map 

Z = ie (2 ^ )w (16) 


For a primary held with h k = h k (of the form cp k (w,w ) = (p k (w)(p k (w)), this procedure gives 
0 (for other primary helds, the one-point function vanishes) 


f dz \ f Qz \ 

(0*(o-,t)) dyn = {<f>k(w,w))gtr = GW 1 TT— {<f>k(z,z))uHP 


= a k (e 2 ^ + e 


\dw J \dw ) 

~ 2 ^)~ 2hk ~ a.e-^ + ..., 7 < 0) = 2vrA k /p = 4nh k /(3 (17) 


1 The subscripts str, cyl will denote a ‘strip’ and a ‘cylinder’, respectively. 
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We have used the following result for the one-point function on the UHP: 

(M z , z ))uhp = A k ((pk(z)(fl(z'))uHP = A k (z - z') 2hk , h k = h k , z'— z (18) 

which follows by using the method of images where the antiholomorphic factor of (j) k (z,z) 
on the upper half plane at the point (z, z) is mapped (up to a constant) to the holomorphic 
ip* k El on the lower half plane at the image point (z', z!) with z' = z,z' = z [23], 22]. In the 
above a k ,A k are known numerical constants. Note that 

Z = ie 2n A+ ir )/P = ie M°-t )//3 ; z ' =~ = _ ie 2n(a-ir)/f = _ ie 2n(v+t)/f) ( 19 ) 

so that in the large time limit we have 

t —> oo => z —> 0,z —» —ioo. (20) 

The second, thermal, expectation value in (fl5]h for ji n = 0, corresponds to a cylindrical 
geometry in the ru-plane, with r = 0 identified with t — /3. By using the same conformal 
map (fTUj) this can be transformed to a one-point function on the plane. For a primary held 
the latter vanishes. Hence (0 is trivially satisfied. 

For a quasiprimary held cf > k , its conformal transformation generates additional terms, 
including possibly a c-number term c k (e.g. the Schwarzian derivative term for <p k = T ww ) 
and generically lower order operators. The c-number term does not distinguish between a 
plane and an UHP. This leads to the following overall result (for fi n = 0): 

(0/c (^"))eqm ( < / > fc(^) ^))cj/Z Cfc, 

(0 fc (a, t)) dyn = (< fa{w , w))str = Cfc + a + ..., yf ] = 27rA fc //3, (21) 

where A*, now is the scaling dimension of the minimum-dimension operator in a T(zi)4>k(z) 
OPE. This is clearly of the general form (J6]) for = 0. 


We now turn to a discussion of these expectation values (1131) in the presence of chemical 
potentials /i n ,n = 3,4,..., as in ([5]) and (j2J). We will denote the new conserved currents as 
W n (w ) and VV n (tc), n — 3,4,.... The conserved charge, W n , is dehned as 


W n = 


2t r 


W TT T da = 


2vr 


(I'dwi W„ ! V' 1 j + (~i) n dw 1 W n (w 1 )) 


( 22 ) 


Here the contour T is taken to be a r = constant line along which dw\ = dw 1 = da. Under 
the conformal transformation (TT6li to the plane/UHP, the holomorphic part of the contour 
integral becomes 


wy„ OI = i-1 


n —1 


2n \ f3 J J Tl 


dzi 


K2J 


u X(«i) + ^2 


«n,n-2m‘l 2 ™ 1 hV n _2m(^l) 


m=l 


(23) 


2 We distinguish from <p k to allow for charge conjugation. 
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where the a n ^ n _ 2 m denote the mixing of W n (z i) with lower order W-currents under conformal 
transformations [21]. The contour r x is an image of the contour T onto the plane. The 
expression for the antiholomorphic part W n \ an tihoi is similar. 

As mentioned before, in this paper we will regard the W n as conserved charges of a W- 
algebra, although the results we derive will be equally valid as long as these charges, together 
with H , form a mutually commuting set, and the currents {YJ n {w), VJ n {w)) are quasiprimary 
fields. 


2.1 One-point function on the cylinder with chemical potentials 

For simplicity we first consider the equilibrium expectation value in (fT5]h Unfortunately, 
unlike the thermal factor above, the factor e~ > lnWn in ([3]) cannot be dealt with in terms 
of a conformal map. We will, therefore, treat this factor as an operator insertion, and write 


{M°))e qm = Tr (0*(w, w)p eqm { l 8, fin)) = - ^^ > W ))cyl ( 24 ) 


We will now illustrate how to compute this for a single chemical potential, say /x 3 , using 
perturbation theory Feynman diagrams!* 3 ! 


lA 


(M^,w)r cyl = (M^^)) cyl - Ps(WsMw,w))^ n + f ( WtWsMv,*))™ + Ofa 


(25) 


The first term in the above expression is the constant Ck that we already encountered in (1211) . 
For a holomorphic primary held cf>k , the second, 0(/i 3 ), term, transformed on to the plane, 
gives 


[ dzi + HT [ dz i zKmtzjMz))^ 

Jr! Jr i 

(26) 

Here we have used the contour representations (l22j ) and (123|) . The correlator inside the 
second integral obviously vanishes (it factorizes into a holomorphic and an antiholomorphic 
one-point, functions, leading to a vanishing connected part). The first integral vanishes unless 
4>k = W 3 (this uses the orthogonality of the basis of quasiprimary fields). In the latter case, 
using 

(WM)W 3 {z))c = 

the integral evaluates to c/(9CU 3 ); combining with the factor of z 3 outside (hk = 3 in this 
case) we get a ^-independent constant, as we must, because of translational invariance on 
the plane. With an antiholomorphic primary held 0*., the calculation is isomorphic. For a 
primary held with nonvanishing hk, hk the result vanishes. For quasiprimary 0*., as well as 
for other W n charges, the conformal transformation to the plane additionally generates lower 


2vr 


<W 3 0 fe H>-r = 


hk 


13 The superscript conn denotes ‘connected’. 







order operators (see, e.g. (125]) )). each of which can be dealt with as in (125]) . The result is a 
finite constant which we will denote as 


(W n (p k (w,w)) = c n , fc 


(this will be non-vanishing only for special choices of 0*,, e.g. 0^ = W n ). As explained above, 
for n = 3 and 4>k(w,w) = W 3 (w), c n ^ = — 27tc/(90/3 2 ). 

In a similar fashion, the 0(/i§) term in (125]) can be transformed to the plane. Again, we 
present the explicit, expression for the simple case of a holomorphic primary held 0*,. 


= (||) V* 


i b I dzi i dz 2 {W 3 (z 1 )yV 3 (z 2 )(j)k(z)} c f nn z\z\ 


+ (-if [dz\ [dz 2 (W 3 (z 1 )W 3 (z 2 )(j)k(z)) c f nn i\zf + / dz x I dz 2 {W 3 (z 1 )W 3 (z 2 )(fik(z)) c f nn z\i^ 


'r 2 

V conn —2 —2 


' Ti 7T 2 


' Ti J r 2 


/ d£ 2 (W 3 (zi)W 3 (z 2 )0 fc (;0)c onn ^ 


'Tj Jr 2 


(27) 


For holomorphic quasiprimary 0*., additional, similar, terms appear due to the generation of 
lower order operators under conformal transformation to the plane. Only the holomorphic 
correlator survives (as in the 0(fi 3 ) calculation). Thus, e.g. if 0*, = T(z ), the stress tensor, 
we have 


(m(zi)W 3 (z 2 )T(z)) c = 


c 

(o - z 2 f(z 1 - z) 2 (z 2 


z) 


2 


Again, after performing the integration over Z\ and z 2l we obtain a ^-independent constant, 
as we must. The analysis of more general fields 0& and two arbitrary IF-charges is straight¬ 
forwardly generalizable. The result is a finite constant (can be zero for a particular 0*,) which 
we denote as 

(hF m hF ra 0fc (w, re)) c mn y. 


Note that in (1271) the result does not depend on the location of the contours T 1 ,T 2 on the 
plane, since the W -currents are conserved. 

Summarizing, we get 


(0fc(re, re)) c ^ Ck f ^ Jl n C n k T j ^ ^ i^m^n Cmn,k T 0(n ) 


m,n 


(28) 


2.2 One-point function on the strip with chemical potentials 

Similarly to the previous subsection, we will treat the /i-deformations in ([5]) as operator 
insertions: 


{M a ^))dyn = (0(t)|0fc(a)|0(f)) 


(e Y.nVnWn/4: 0 fc ( W;W ) e E„ W n / A ^ 

( e -EnMnW„/2) str 


(< h(w,w)Y str 
(29) 
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As before, we begin by illustrating the calculation of this quantity with the simplest case of 
a single chemical potential /r 3 , using perturbation theory Feynman diagrams: 

((j) k (w,w))f tr = {<f> k (w,w)) str - ^({ W 3 ,(f) k (w,w)})^ n 

+ Of) Mw, u.)})r” + 2(1 V 3 Mw, w )W 3 )Z n ) + O(pl) ( 30 ) 

The { , } denotes an anticommutator. The operator ordering implies the following: when W 3 
appears on the left of cp k (w,w), e.g., in (W 3 (p k (w,w)), the integration contour (122j) for W 3 
on the strip lies above the point ( w , w); similarly when W 3 appears on the right of w ), 
e.g. in {(p k (w, w)W 3 ), the contour for W 3 is below the point ( w,w ). 

The first, /i-independent, term in the above expansion is already calculated in (j2Tj) . 

2.2.1 0(fi n ) Calculation 

As before, we find it convenient to use the conformal transformation (fTGlh The correlator on 
the strip then reduces to that on the UHP, as in the ji = 0 case before. For a holomorphic 
primary field <p k , this gives 

z 3 [ dz ! + H ) 3 [ dz x zKmiz^M^VHP 

■/iy Jr 1 

(31) 

where the operator ordering explained above implies that the contour Ti lies to the left of 
the point ( z,z ) on the UHP. Now, in the analogous calculation (126|) . the second connected 
correlator on the complex plane vanished because of factorization into one-point functions. 
Correlators on the UHP are, however, related to those on the plane by the method of images 
(an example of which we saw in (1T8]P . In particular, VV 3 at the point (z\,z\) on the UHP 
becomes the holomorphic operator W 3 = — W 3 on the LHP at the point (z(, z[) with z[ = z 3 
[23U22]. The contour Ti gets mapped to its mirror image r( on the lower half plane. With 
this, we get 


2t r 


^ hk 


p 2 


2tt 


(w 3 M^))Ttr n = ^z hk 


P 2 ' 


Hi+r' 


dzi «?(W 3 (^i)0fc(2))g” 


(32) 


On the complex plane, the contour V 1 on the UHP can be deformed to r( on the LHP, hence 
the two contours simply yield a factor of 2. In fact, combining with the other ordering, and 
applying a similar reasoning, we get an overall factor of 4. Thus, combining with results 
from Section I2TT1 we get, for holomorphic primary fields 

(33) 

A similar statement is true for an antiholomorphic primary field. 

Let us turn now to primary fields (j) k (w,w ) with h k ,h k 7 ^ 0 (of the form (j) k (w,w ) = 
tp k (ui)tp k (w), as discussed before in the context of (fT5|P . In the cylinder calculation in 
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Section 12.11 the /i- correct ions for these vanished. In the present case, they are non-zero for 
operators of the form (p k (w,w) = <pk(w)<p k (w), with h k = h k (as in (jT7|) ). After conformally 
transforming to the UHP, we regard (p k on the UHP as (p k at the image point on the LHP 
(up to a constant). Combining with the arguments used for the holomorphic operators, we 
eventually get 


{{W 3 ,(j) k {w,w)}) 


conn 

str 


str 

h M= [ 


?7T 

= z 3 -(zz) h I 3 (z,z'), 


ri+ri+fi+f; 


(34) 



Figure 1: Various contours needed to compute the W n insertions in (13U1) . At late times, the insertion 
of each contour, irrespective of the position of the contour, amounts to insertion of a given factor linear 
in t. This allows to resum arbitrary orders of arbitrary WL-charge insertions, leading to the exponential 
time-dependence as in ©. See Figure [2] for more. 

The ratio of correlators inside the integral is given by 


(W,( Zl ) Vt ( z )^( 2 ')>r7(^(bvl(^')>r" = 937 —(35) 

where q 3 is the W 3 -charge of the field cp k . Integrals of the kind (1341) are discussed in detail in 
Appendix IA.21 The final result (see f|7ip ) is that the 0(/i 3 ) correction, in the long time limit 
([HD , is given by (using that all four contours Ti, Ti, r(, T 7 , contribute equally, cancelling the 
1/4 in -/i 3 /4) 


(0fc(o-,i))dyn = a k e 2,rA *A/3 ^1 - Q 3 , k n 3 (— + constant) + 0(/i 3 )J + ..., 

Q 3tk = i 3 2q 3tk (27r), fl 3 = p, A k = 2 h k (36) 

Up to 0(/i 3 ), it agrees with (El). 

In case of a quasiprimary field 4> k (w,w), it mixes, under conformal transformation to the 
plane, with lower dimension operators. The most relevant operator among these, which is of 
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the form tp k (z)(p k (z), is then to be used in (j5jj) for obtaining the dominant time-dependence; 
in that case A*,, Q^ k refer to this operator (rather than to the original 0*.). 

For higher W n charges, the currents W n (w) are typically quasiprimary, and hence they 
mix with lower order W m (z) under conformal transformation to the UHP. Thus the 0(p n ) 
correction to the dynamical one-point function (c t>k)dyn is a linear combination of terms of 
the form (jfigj) (weighted by a set of coefficients a n>m , as in (1H71) below). Collecting all this, 
the O(n) correction with all chemical potentials is given by 

( 0 fc (* 2 ") t) ) dyn 

Hn 
Qn,k 


Note that for W 3 deformations, the expression for Q 3 as in (j5S|) corresponds only to the first 
term in the above series expression for Q n . This is because the W 3 current is a primary held 
and does not mix with any lower W current under a conformal transformation. From n = 4 
onwards, the additional terms in Q n ,k s represent the mixing of W„ currents with W n -2m 
under conformal transformations. 


= a k e~ 2nAkt/ P 




f 1 _ + constant) + 0 (/i 2 ) j + ..., 


— Rn _i 1 — h k + h k — 2 h k 


f3 n 

Ln/2—1J 


2 ^ a n , m i n - 2m (27r) n - 2m - 2 q n - 2 m,k 

= i n ( 2 ir) n - 2 2 q n , k + i n ~ 2 ( 2 Tr) n -*a n , 2 2 q n _ 2 , k + ..., 


m =0 


(37) 


2.2.2 Higher order p-corrections 

Let us hrst consider that 0(/x 2 ) correction: 

(Mu.'wKrr = {jh ^r-(({W n W n ,Mw,w)}) c ,r + 2(W n Mw,w)W n )™ n ) (38) 

Again, for holomorphic (or antiholomorphic) primary fields (j) k (w), it is straightforward to 
generalize (133|) to this order. 


= %r(W„W n Mw))% 


2 ! 


(39) 


For a primary held of the form <f) k (w,w) = cp k (w)ip k (w), proceeding as in the previous 
subsection, we get 


1 { 277 \ ^ 

(4>k{w))'i™ n I 2 " - [Qn.k^nt — ) + At); (constant X t + constant) + 


(40) 


The essential ingredient in this calculation is 


I nm (jH,z'\T 1 ,r 2 ) = / dzi Z " 1 / (lz 2 Z™ 1 f nm (z U Z 2 , Z, z'), 
dry dr 2 

(Wn{z 1 )W m (z 2 )< Pk (z)<Pi(z')) C C° nn 


fnm(z 1 ,Z 2 ,Z,z') = 


{Mz)<P* k (z')) c c 


(41) 
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By repeating the strategy of (1751) . we get 


Coefficient of [log (-s') - log(-z)] 2 in I nm (z, z'\T 1 ,T 2 ) 


= Residue. 1=z 


Residue. 2=2 


(W n (zi)W m {z 2 )ip k (z)ipl(z')) c ^ 

(Mz)vl(z')) c r n 


Qn,k Qm,k 


(42) 


where we have first used the Wmiz 2 )(pk(z) OPE, and then the W n (zi)<^fc(z) OPE. In a manner 
similar to that in Appendix IA.21 we conclude the following structure of I nm (z,z'): 

Inm(z , PRi, r 2 ) = q n> k q m ,k{{ log(-z') - log(-z)] + constant) x ([log(-z') - log(-z)] + constant) (43) 

Note that at late times f > ,3, ([log(— z') — log(— z)\ —> 2(2nt)//3 and dominates over the 
constant term (the precise sense is that of (ITS]) ). Similar to Appendix IA.21 the 4 x 4 = 16 
locations of the contour-pairs (Pi, r 2 ),(Pi, r' 2 ), (Pi, r 2 ),(ri, H,), ...., all contribute equally, 
therefore converting (/i„/4)(/i m /4) —$■ /i n p m . Combining all these, we get (ITU]) . The charges 

S that are dehned by the W n <p OPE and appear in (1421) . get multiplied by some constants 
and shifted by lower W n - 2 k charges to give the Q n in (1T7D . as in (11771) . 


Arbitrary orders and Exponentiation: 

It is straightforward to generalize the above 0(ft, 2 ) calculation to higher orders in the per¬ 
turbation in chemical potentials. Thus, at the order H*=i /Pi,, there are r insertions of 
W-currents, leading to integrals of the form 


Inin 2 ...n r % | Tl j T2) TV) — / dZ\ Z^ I dZ 2 Z. 


ni-1 I ^n 2 -l 

ri J r 2 


dz 2 z? r 1 fn 1 n 2 ...nr(z 1 ,Z 2 , Z r )Z,z'), 


t „ _ (W ni (z 1 )W n2 (z 2 )...W nr (z r )g} k (z)g}l(z'))^ 

Jniri2...n r •••? Z ) 


(<Pk( z )'Pk(z')) C c 


(44) 


Again, repeating the strategy of ([75]) . we get the following leading (viz. (log) r ) contribution 
(see (ITS]) for the definition of the leading-log contribution) 

Coefficient of [log(-z') - log(-z)] r in J Ttl „ 2 ...„ r (z, z'\T 1 , r 2 ,..., T r ) 

( (W ni (z X )...Wn T _ x (z r - 1 )>V nr ( 2 r ) W (z)^(«'))g m 


= Residue. 


...Residue, r _ 1=z < Residue Zr= 




(v k (z)v* k (z')) c‘ 


— Qni,k”‘Qn r —i,k Qn r ,k 


(45) 


where we have hrst used the W nr (z r )(pk(z) OPE, then y^ nr -A z r-'i) ( Pk{ z ) OPE, etc. As in 
the 0(/i 2 ) calculation above, we obtain the following behaviour at late times 

In 1 n2-..n r (z,Z | Tl. T 2 , ..., T r ) 

= 9ni,fc-9n r _i,fc Qn r ,k ([log(-z') - log(-z)] + constant) x ... x ([log(-z') - log(-z)] + constant) (46) 

N --v---' 

r terms 

The two equations above show that the leading log contribution to (21 from, every contour 
integral of the W ni current contributes the factor q n . [log(— z') — log(—z)]. This is the first 

14 Each W n current comes with a factor of ^ ^ j , as in (l23l) . 
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basic ingredient for the exponentiation we are going to find. Furthermore, it is easy to see 
that the leading log contribution is the same irrespective of where each contour Fj is placed 
(out of 4 possible choices, e.g. Ti, T^, r i, in Figured]). As before we must combine the 
contribution of all positions of the contours, which, therefore, amounts to multiplying the 
result for (144|) by 4 r which converts the original coefficients coming from exp[—/i n PF n /4] 
as follows 

nk ^/ 4 nr=i ^ 

v\ r\ 

This is the second basic factor leading to the exponentiation. Combining all these, and 
incorporating some additional constants (see footnote [TTj) we get the following, leading, 
order (fJ, ni ...fjL nr ) contribution 


(M w )) 


conn | 
str I 


| r 



0(ii r t r ~ l ) 


(47) 


Once again, the constants Q n are related to the q n as in 113?]) ) in a manner similar to the 
O(fi) and the Ofifi 2 ) calculation above. We note that the leading log contribution used in 
this paper can be isolated by considering a scaling 

—* 0, t = — —> oo, such that fi n t = constant. (48) 

r 

The second term in (147p . or for that matter, in (HU]) . is subleading at large times in the sense 
of this scaling. 

Using the above results, we now have, for primary fields of the form (f) k (w, w) = <p k (w)(p k (w) 


(<f) k (w,w)) str = a k e p 


2nt 


1 - > , h'n Qn,fe(-y- + const) 

1 v~~v „ _ , 27 t t . .27rt 

+FT / .hnUm Qn,k{— + const) Q m ,k{— + const) + .. 


2 ! 

1 

t\ 


p 


} ' TT hriiQn ir 


{m}i=i 
ake 


/ 

\ 



27 it 2iit 


( 

——b const)•••(——b const) 

jj (J 

+ ... 

V 

V 

r terms / 



— aup -2nt/P( A k+T,„i i nQn,k+0(p, 2 )) _ a e ~7kt 


(49) 


where we have absorbed some constant factors in a k . 7 k is given by (JT]) ; Q n ^ k are the shifted 
W n charges of as dehned in fl37j) . The proof of the above equation for general quasiprimary 
operators fi k (w,w) works out much the same way as in case of the 0 (/x) terms, as discussed 
in Section 12 . 2.11 We emphasize that it is only the leading contributions at large times which 
we have proved here to exponentiate. Thus, we do not claim that the constant terms marked 
“const” in the above equation are all the same. As we have remarked before, the leading 
contributions can be isolated using the scaling mentioned in (1481) . 

The schematics of the above calculation is explained in the Figure [21 
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Figure 2: The schematics of the calculation of the one-point function. The first term represents the zero- 
order boundary Green’s function (1171) without chemical potentials (the shading indicates the boundary of 
the upper half plane). The second term represents the 0(p, n ) correction, which involves one insertion of a 
Wn-charge (which is an integral over the zi-contour. As explained in the text, at long times, this insertion 
amounts to multiplying the zero order term by a term of the form / n log(;z), where /„ is described in (1371) . 
The third term represents insertion of two such W-charges; as we explained in the text (see m and below), 
each insertion again amounts to multiplying by the factor mentioned above, along with a factor of ^. The 

pattern continues, to ensure an exponentiation to Gq(z) , as in (l49l) . Since at long times G o(z) ~ e -7 * 

(see (ED), and z ~ e 27rt /' s , adding the chemical potentials amount to a shift of the exponent —> 7 k as 

in ©. 


3 Calculation of I(t) 

Let us rewrite the expression for the thermalization function I(t) (jSJ) in the form 

I (t) = Z sc / \J Z SS Z CC = Z sc j yj Z SS Z CC , 

Zsc = Tl i.Pdyn,A y^)Peqm,A (/^, h) ), Z sc Z sc / (^Z S Z C ^ 

Zss = Tf (pdyn,A(t)Pdyn,A (t)), Z ss Z ss /Z s , 

Zqc = Tl (Peqm,A (/3, P)Peqm,A^Pi P) ), Z cc Z cc /Z c , 

Z s = Tr (p dyn (t)) = (ipolipo), z c = Tr (pp^) (50) 

In Appendix [B] we explain how to compute I(t) using the short interval expansion, valid 
when the length of the interval / is small compared with the other time scales /3 and t in the 
problem. We reproduce the main formula (1751) for our purpose, where we explicitly denote 
the dependencies on the length l of the interval, the inverse temperature /3 and the chemical 
potentials /i (the dependence on f3 on the RHS is implicit; the one-point functions depend 
on both /? and p — see Section [2]). 

Zse(l,P,fj) = C 'fc 1 »fca(0(‘^fci( w l» ^>l))^T-(^fca( w 2, ^> 2 ))^, 

ki,k 2 

Zss(l,P,V) = C'fcl M (0 (• Wl , Wi ) ) ^tr (0fc2 (' w 2 , W 2 ) ) sir, 

k\,k 2 

Z cc {l,P,n) = G fcl,fc2(O(0^(^l,^l))^(0fc 2 (w2,W2))^Z (51) 

ki,k 2 

It is understood, for the logic of the short interval expansion to go through, that all contours 
which represent insertion of the W-charges (see Fig[T]) are drawn outside of the small disc-like 
region of both sheets of Fig QJ 
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3.1 Proof of thermalization 

Using the short-interval expansion above, and the long time behaviour of one-point functions 
from Section [2]), it is easy to prove that the system thermalizes in the sense of (TTOl) or CO- 
To prove this, note that it is only the holomorphic (or antiholomorphic) fields (pk which 
possibly have non-zero expectation values in the long time limit (l20|h For these fields, the 
one-point functions on the cylinder and on the strip agree (see (OTT) . (1331) . (IHH ). By virtue 
of (T5T]1 , we therefore have in the long time limit Z sc = Z ss = Z cc . Hence using the expression 
(l50|i for the thermalization function we get I(t —> oo) = 1 which proves (HOjl and consequently 
00 - 

The above-mentioned equality of one-point functions between the strip and cylinder ge¬ 
ometries for holomorphic (or antiholomorphic) fields imply the same for the conserved W n - 
(or VV n )- currents. This, therefore, proves that 

(m\Wn\m) =^iW nPeq m) (52) 

Note that in proving this, we have used the correspondence ([3]) between the parameters of 
the initial state and the putative equilibrium state. The above equation, therefore, proves 
the correspondence ([4]). 

3.2 Thermalization rate 

To evaluate the rate of approach of I[t ) to its asymptotic value 1, we organize the terms in 
Z sc , Z ss , Z cc as follows 

z,c = c 0 ,o( 1 + sn, sr = + <0„}y + Y.C.MXMI, 

a ab 

Z„ = C„,„(1 + 5" + s 2 "), 5" = 2 V q,„(0 a >; tr + V sr 

a ab 

= c„,„(i + sn. sr = 2 £ c aS {<t,x„+Y. 

a ab 

where a,b ,... denote descendents of the identity operator, k labels other primaries (than the 
identity) and their descendents. C = C/C 0)0 . 


k 

(53) 
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3 . 2.1 /i = 0 

Let us first consider the case of zero chemical potentials. Using the results in Sections [21 and 
Appendices [A] and IB.11 we get 


S a 2 ‘ = 


<xt — 


—Gy 1 / 2 

(l + o(if ) 

+ C 


+ 0 (e~ Sni 

—(3(p/“ 

(l + 0(7) 

+ 2 d r j' r f'l c 

~ 87ri (1 + O(0 2 ) 

+ 0 (e~ 8ni ) 

^2 \a k l Ahk e~^ hkt/p (l + 0 (1) 

k 

2 ) + 0 (e~ 12nhki ) 


— cltI 2 

(i + o(p=) 




C7T 2 

A T 7T 4 

a\ t 

7T\ ih k 


24 ' “ TT “ 8c 

a k = — 
n k \ 

2 / 



(54) 


To this order, it is easy to see that the contribution to I(t) from descendents of identity, 
demarcated by cqr, a T f, vanishes. The leading contribution to I(t), demarcated by a*,, occurs 
only in Z ss and comes from ((4> m (z, z)) str ) 2 for which h k is the minimum (= h m ) (this could 
be a held which appears after a conformal transformation of the original quasiprimary held). 
The time-dependence shown of S™ comes from (fTTjh Using this, we get 

I{t) = 1 - aexp[-2^t] + ..., yj® = 2vrA m /(3 (55) 


This is of the form ()9|) for fi — 0. with 

0=^(1)“" (56) 

The discarded terms in (1^51) are faster transients. This proves (ED for zero chemical potential. 
This result has already appeared in 


3 . 2.2 // ± 0 

The generalization of the above result to the case of non-zero chemical potentials is straight¬ 
forward. Once again, the dominant time-dependence arises from ((</> m (z, z))^ tr ) 2 in the 
or Z ss . The time-dependence (ED follows by using (BED in ■ 


3.3 Properties of Q 

From the asymptotic behaviour (ED of the thermalization function we indicated the asymp¬ 
totic behaviour (fT3l) of the dynamical reduced density matrix Pdynif)- By using the long 
time behaviour of the one-point functions (ED, we can easily deduce the following dominant 
behaviour of overlaps of Q with various quasiprimary fields at late times 

Tr (Q(p k (t)) oc Tr (Q<f> m (t)) —> constant. 

15 0ur exponent differs from Cardy’s value by a factor of 2. 
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4 Decay of perturbations of a thermal state 


We found in the previous sections that the long time behaviour of the reduced density 
matrix pd yn ,A{t) resembles that of a thermal ensemble plus a small deformation which decays 
exponentially. We will find in the next section that the thermal ensemble (or more accurately 
the generalized Gibbs ensemble) corresponds to a (higher spin) black hole geometry in the 
bulk. The small perturbation of the equilibrium ensemble is thus expected to correspond to 
a small deformation of the black hole geometry. Consequently, the exponential decay of the 
deformation in the CFT should correspond to a ‘ringing-down’ or a quasinormal mode in 
the bulk. 

We will address the above issue in the next section which deals with bulk geometry. 
However, in order to make the correspondence of the above paragraph more precise, in this 
section we will directly present a CFT computation of the decay of a perturbation to a 
thermal state. Note that this computation is, in principle, different from the exponential 
decay of the one-point function in the quenched state, (J6]). However, what we will find is 
that the long time behaviour (GU) of an operator 0 fc (O ,t) in the quenched state is the same 
as that of its two-point function (1571) in the thermal state (j3J) (with chemical potentials). 
The latter measures the thermal decay of a perturbation and is more directly related to a 
black hole quasinormal mode. Throughout this section, we will assume that the conformal 
dimensions of </>& satisfy hk — hk- 

We define the thermal two-point, function as @ 

G+(t, 0 ;/ 3,/4 = |Tr(0 fc (O,t)0 fe (O,O)e-^-^^) (57) 

By the techniques developed in the earlier sections, a computation of this quantity amounts 
to calculating the following correlator on the plane 

{(,t> k {z,z)(p k (y,y)e^^^ Wn ), z = ie~ 2nt/fi ,z =-ie 2 ^ t/f3 ,y = i,y =-i (58) 


where the p n -deformations are understood as an infinite series of contours as in the previous 
section. 

For p — 0, the above two-point function is standard. Including the Jacobian of transfor¬ 
mation, we get 

G + (t, 0; /3, 0) = (j) 4h * [(ie- 2 ^ -i)(-ie 2 ^ + i)]~ 2hk ^ const e ~ 2 * tA ^, (59) 


which clearly matches the long time behaviour of the one-point function (JSJ) in the quenched 
state for p = 0. Here A*, = 2 hk- 

I 11 the above, we considered the thermal Green’s function for two points which are both 
at the same spatial point a = 0. It is easy to compute the Green’s function when the two 
points are spatially separated by a distance Z, say with <xi = Z and a 2 = 0. We get 


G+(t, Z; /3, 0) = -Tr(0 fc (Z,Z)0 fc (O,O)e-^) = 


2vr 

IT 1 


,irl/0 


1 4/ifc 


f const e 27rtA fc/^ (t — Z) (3 
* < const e - 2 ^ A GJ, (Z - t) » P 


^ ie 2n(l-i)/P _ ^_ ie 2n(l+t)/P + 2h * 

( 60 ) 


16 We use the same notations as in 
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The coordinates of the two points, in the notation of (158|) are modified here to z — ie 2n ( l ~ t ' > /P , 5 = 
_i e 2 Tr( l +t)/p^ y — i^y — —I. The prefactor with the square bracket comes from the Jacobian 
of the transformation from the cylinder to the plane. The behaviour of the Green’s function 
is shown in Figure |3j It is important to note that the exponential decay, found in ((6]) shows 
up only for time scales t l. 



Figure 3: Plots of the thermal Green’s function G+(t,l;j3, 0) for /3 = 27r, = 1.5. The curve on the left 
(blue) is for l = 6 , and the curve on the right (orange) is for l = 8. Note that the exponential decay in time 
occurs for times larger than l. 

The effect of turning on the chemical potentials can be dealt with as in the previous sec¬ 
tions. At 0(/i n ), we will have, as before, a holomorphic contribution and an antiholomorphic 
contribution. The former is proportional to 

{(t>k{z)(j)k{y)) x J dziz™- 1 (W n (zi)(j)k(z)(j)k(y)) (61) 

As we see, the structure of the integral is the same as in the previous section. As before, 
logarithmic terms appear in the above integrals which give the leading, linear, f-dependence. 
Similar remarks also apply to the antiholomorphic contour. Since the calculations are very 
similar to those in the previous two sections, we do not provide all details. By resumming 
the series over the infinite number of contours, we find in a straightforward fashion that 

G+(t, 0; n) G+( 0, 0; /J, 0)6(/i)e“ 7fct (62) 

where b(fi) is time-independent, and is of the form 6(/i) = 1 + 0(/i). This long time decay is 
the same as that of the one-point function (J5J) in the quenched state, as claimed above. For 
points separated by a distance Z, the above exponential decay shows up for t 3> Z, as in (l60]h 
In the above, we have discussed the two-point, function in real space. It is straightforward 
to convert the result (1601) without chemical potentials to Fourier space, which develops poles 
at 

i 2 tt / » \ / \ 

| M=0 = k + 2m), m = 0,1, 2,... (63) 
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Our results in fl5} can be interpreted as a shift, caused by the presence of the chemical 
potentials of the dominant pole o| M=0 t° 

2 n _ 

0 ^ q k “b ^ ^ nQn,k ) (64) 

" n 

where the notation is the same as that of (J6]) . In this paper we will not address the question 
of the shift of the subdominant poles u)k,m (for m = 1, 2,...) due to chemical potentials (the 
current status of these can be found in [26] [9, HU]). 

Two-point functions of the kind fl-57]) . for a single chemical potential /i 3 , and up to order 
/ig, have appeared earlier in p) (calculations up to 0(/i|) have appeared in PH])- What we 
find in our paper is that at large times, the perturbation series in /d n , up to all orders in all 
chemical potentials, can be resummed, to yield the leading correction to the thermalization 
rate in the presence of chemical potentials. 

At a technical level, the one-point function in the quenched state corresponds to a one- 
point function in a geometry with a boundary, and for operators considered here, these turn 
into a two-point function on the plane, by virtue of the method of images. The thermal decay 
naturally involves a two-point function on the plane E3 and agrees with the above two-point 
function at late times. 


5 Holography and higher spin black holes 

Zero chemical potential: As remarked in the Introduction, a global quantum quench 
described by an initial state of the form (J5]) , for large central charges and zero chemical 
potentials, has been shown in HZl US] to be dual to one half of the eternal BTZ (black 
string) geometry, whose boundary represents an end-of-the-world brane. 

In an independent development, it was found in 123 that the quasinormal mode of a scalar 
field <1 >k(cr,t,z) of mass m in a BTZ background (dual to a CFT operator 4>k of dimension 
Afc = 1 + \/l + m 2 ) is of the form exp[— 2nAt//3] at large times. This time-dependence 
agrees with the CFT exponent in (IUH exactly. This shows that the exponential decay of a 
CFT perturbation to a thermal state corresponds to the decay of the corresponding scalar 
field in the bulk geometry. This result has been extended to higher spin fields in the BTZ 
background in [28] . 

Non-zero chemical potentials: In case the CFT has additional conserved charges, in 
particular if it has a representation of a algebra (and consequently the hs(A) algebra 
US]), then the bulk dual corresponding to those conserved charges have been conjectured to 
be the conserved higher spin charges of higher spin gravity. In particular, [2U] have shown 
that if one interprets the grand canonical ensemble (J3]) (more generally, the GGE) in the 
framework of an hs(A) representation, then the bulk dual corresponds to a higher spin black 
hole. 

17 Actually the thermal calculation involves a product of two such factors, one holomorphic and the other 
antiholomorphic, but one of the factors just gives an overall constant and only one factor leads to the 
important time-dependence. 
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Thus, we would like to conjecture that the bulk dual of the quantum quench with chemical 
potentials, would correspond to a gravitational collapse to a higher spin black hole. 

As an important consistency check, by analogy with the case with zero potential, in the 
present case too, the leading quasinormal mode (QNM) of a scalar field <f>fc(cx, t,z) should 
have a time-dependence given by (1621) . Following the results in |26j (see also p, uni ei]) Efl 
we find that at late times t f3 the QNM for the hs(A) scalar field $ + behaves, up to 0(/i3), 
as e~ lulk ’ ot , where 


Wk, o 



1 + A + F3g(l + A)(2 + A) 


(65) 


where the index k here refers to the operator 0*. dual to the scalar field <3> + . Noting that for 
this operator we have A*, = 1 + A, and Q^k = |(1 + A)(2 + A) [9, 10], we see that the QNM 
frequency cok,o agrees, to the relevant order, with the pole (1641) of the thermal 2-point function 
which, in turn, is related to the thermalization exponent by the relation cJk,o = with 

7 /c given in (0). 


6 Discussion 

In this paper, 2D conformal field theories were considered with additional conserved charges 
besides the energy. We probed non-equilibrium physics starting from global quenches de¬ 
scribed by conformal boundary states modified by multiple UV cut-off parameters 0. It 
was found that local observables in such a state thermalize to an equilibrium described by 
a grand canonical ensemble (J4]) with temperature and chemical potentials related to the 
cut-off parameters. We computed the thermalization rate for various observables, including 
the reduced density matrix for an interval. It was found that the same rate appears also in 
the long time decay of two-point functions in equilibrium (see (0) and ffl4|i b In the context 
where the number of conserved charges is infinite, and they are identified with commuting 
Woo charges, the equilibrium ensemble (a generalized Gibbs ensemble, GGE) corresponds to 
a higher spin black hole [20] • We found that the thermalization rate found above agrees with 
the leading quasinormal frequency of the higher spin black hole; this constitutes an addi¬ 
tional, dynamical, evidence for the holographic correspondence between the global quenches 
in this paper and the evolution into the higher spin black hole. 

One of the main technical advances made in this paper is the resummation of leading- 
log terms at large times, presented in Section 12.2.21 which leads to exponentiation of the 
perturbation series, leading to the thermalization rate, presented in (0, (H9]h as a function 
of chemical potentials. This allows us to also compute the effect of chemical potentials on 
the relaxation times of thermal Green’s functions. Another technical advance consists of the 
computation of the long-time reduced density matrix 0, using a short-interval expansion, 
which allows us to prove thermalization of an arbitrary string of local observables. 

One might wonder whether the results presented in this paper are tied to the use of trans- 
lationally invariant quenched states such as 0, whose energy density and various charge 

18 We wish to thank Alejandro Cabo-Bizet and Viktor Giraldo-Rivera for informing us that the difference 
between equation (1651) above and the corresponding equation (4.2) in a previous version of their paper [26] 
was due to a typo, which has now been corrected in the new version of their paper. 


21 





densities are uniform. We will address the question of inhomogeneous quench in a forthcom¬ 
ing paper pH], both in the CFT and in the holographic dual, using the methods of [22| where 
we create an inhomogeneous energy density by applying conformal transformations. It turns 
out ra that if the initial state has inhomogeneities in a compact domain and has uniform 
energy densities outside, local observables again thermalize asymptotically with exponents 
governed by the uniform densities. Other important issues involve local quenches (see, e.g. 
[30]), and compact spatial dimensions. The issue of thermalization when space is compact 
is quite subtle. It has been shown in [TO] that at large times one can have the phenomenon 
of revival (observables effectively returning to their initial values). The dynamical entangle¬ 
ment entropy for a quantum quench in a space with boundaries is an interesting, related, 
issue; we hope to come back to this in a forthcoming publication [31] . 

Acknowledgement 

We would like to thank Pallab Basu, Justin David, Deepak Dhar, Oleg Evnin, Rajesh 
Gopakumar, Shiraz Minwalla, Pranjal Nayak, Arunabha Saha and Tomonori Ugajin for dis¬ 
cussions and Somyadip Tliakur for discussions and collaboration in the forthcoming paper 
[21] . partial results from which are presented here. We would also like to thank Juan Pedraza 
for drawing our attention to Ref. [3] and Alejandro Cabo-Bizet and Viktor Giraldo-Rivera 
for a useful correspondence regarding Ref. m 


A Some details on one-point functions 

Here we collect some additional helpful material on the one-point functions discussed in this 
paper. 

A.l A few explicit one-point functions with zero chemical poten¬ 
tials 

Case k = descendent of identity : In this case, (j)k(w,w) is of the form T, T, or : TT : or 
some descendents thereof. Under a conformal transformation (fTfijh these operators pick up 
a c-number term in addition to a term proportional to the corresponding operator on the 
plane/UHP. We will give some examples to illustrate the calculation 
1. cylinder: In this case 

(T(w)) cy i = ((~^p ~ = -g^2 

(: TT :(w, w)) cy i = <([-^ - ^*T{z)][-^ - ^T(z)]) W 
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2. strip: In this case 


(T(w))str = 
TT :(w,w)) str = 


C7T 

W 2 


CTT 


4n 2 

4n 2 


z 2 T(z ) ))uhp = 


6f3 2 f3 2 


z 2 T(z)}[- 


C7T 


C7T 

47r 2 


= ( T (w))cyl 


6f3 2 /3 2 


Z T(z)] )uHP 


(^) 2 + 4rf (* - ^) 4 - (^) 2 + a TTe Svt/P + ... (67) 


where A T f , a T -f are constants as in (TT7T) and (ITS]) . 

Case k = descendent of other primaries : In this case, 

1. cylinder: The one-point function vanishes as in the case of primaries. 

2. strip: The one-point function can be related to one-point function of primaries which is 
dealt with above. 


A.2 Some details on 0(/i n ) correction to the one-point function 

In this section we will consider the following integrals which arise in connection with 0(fj , n ) 
correction to the one-point function (</>(cr, f)) dyri : 


I n (z,z’\Ti) = / dzi z™ f n (zi,z,z'), g n (zi, z, z!) = / dzi z" 1 f n (z 1 ,z,z l ) 


' Ti 


, , . „ 

Jtl \ 1 : ? / /. _ / . _* / \ norm. ^In,k , 


(<Pk(z)<P%(z')) c 


\zi - z) n (z 1 - z'Y 


( 68 ) 


The second integral on the first line is an indefinite integral. The integrals above can be 
explicitly computed. E.g. 


g 3 (z 1 , z, z') = q 3 ,k[ R 3 ( z , z')(log(zi ~ z ) ~ logfci ~ z ')) ~ 


v/2 


+ 


2{zi - z) 2 2{z\ - z') 2 


z'(2z + z') z(2z' + z) 


(z - z')( Zl - z') (z - z')(z 1 - z) ! 
(z + z'). 


h(z, z '\Ti) = q 3ik [R 3 (z,z')(-\og(-z) + log(-z')) + 3^--(] 


( z ~ z ') > 




(z - z’) 


>\2 


(69) 


Note that I 3 is essentially obtained from the lower limit of the integral, i.e. from —g( 0, z, z'). 
The contour Ti in I 3 specifies which branch of the log is to be taken. In particular 


/ 3 (z,z'|Ti) - I 3 {z,z' |fi) = -27 iiq 3)k R 3 {z,z') 

In the long time limit ()20)) . we get 

I 3 (z,z'\Ti) = I 3 (z,z'\Ti) = 2q 3)k t(2n/f3) + q 3 , k x const + 0(e" 27rt//3 ) 


(70) 


(71) 
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In this equation we have displayed the principal value of the relevant integrals (the discon¬ 
tinuity (1701) tells us the coefficient of the log term or the linear t term). 

However, we would like to understand the above results more simply, by using the 
W n (*i)^fe(*) OPE which is of the form: 


Wn(zi)Vk(z) = q n ,k + 

where tp k ,i{z) is of dimension h k + Using this, we get an expansion for the connected 
3-point function of the form: 


rc-i / \ 

Y a n.i , Lpk,t v„ i + Aguiar terms (72) 

“ Ui — z) n ~ l 

7 — 1 ' ' 


(VV,,faW z )Wth'))c' 

(Mz)‘Pi(z , ))ir n 


Qn,k 


+ 


4 


n, 1 


(z\ — z) n (z\ — z) n x (z — z') 


+ 0(z - z') 


A—2 


(73) 


Performing the integral in 


9n{zi, z, z') = q^k ^logfzi -z]-{n- 1) ^ Z _ ^ 

c 

+ (zi — z+(n— 1)2 log [zi — z\ + ...) + ... 

z z 


The ellipsis in each round bracket represents terms with higher powers of 1/ (zi — z) (up to a 
maximum of (z\ —z)~ n )] successive round brackets themselves are arranged in higher inverse 
powers of z — z'. Using the W n (zi)(pl(z') OPE in a similar fashion and using the symmetry 
property g n (zi, z, z r ) = (— l) n g n (zi, z', z) we can arrive at a general structure 


g n (0,z,z') = q n , k (log[-z] - log [-z'})R n (z, z') + ... 

where R n (z,z t ) = (—1 ) n ~ 1 R n (z t ,z) is of the form P n _i(z, z')/(z — U) n_1 (P n _i(z, z') is a 
homogeneous symmetric polynomial of degree zero). See the explicit form of R n for n = 3 
in (1691) . The omitted terms are all ratios of homogeneous polynomials in (z,z r ) of the same 
degree in the numerator and in the denominator. This implies that we have, in the long time 
limit (1201) 

J n (z,U|ri) = / 3 (z,U|fi) = 2q rL:k {2'K/@)t + q n)k x const + 0(e _27rt//3 ) (74) 


which, of course, agrees with (I7TI) . 

Note that the dominant time-dependence 2q ntk t(2n / j3) comes from the long-time limit 
of the coefficient R n (z,z') of the log terms, which can be read off from the discontinuity 
4(^1^)-4 (z,z'\Ti) (see my Now, the contour f r f dz\ can be deformed to a very 

small circle <fT z dzi around the point z\ therefore the leading long-time behaviour R Z P (z, z') 
can be derived by using the leading OPE singularity in (172]) and computing the residue at 
Z\ = z: 


Coefficient of [log(— z') — log(— z)\ in I n (z,z‘ 

(W n (zi)^ fc (^)^(0'))r 


= Residue 2l=z 


&k{z)ip* k {z')) ^ 


— Qn,kRn ^(-2'j Z ) qn,k 


19 This is the general form; some of the a U) i coefficients may, of course, vanish. 


(75) 
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B Short interval expansion 

In this section we will explain a formalism suitable for computing partition functions of the 
kind that appear in (l50l) . For convenience we will first compute these quantities in Euclidean 
time t = it and later analytically continue back to Lorentzian time. With this, each of the 
expressions Z sc , Z SS1 Z cc is of the form 

Ti(p a, 1 Pa, 2 ) = I Dy?i f D(p 2 8(F[<p 1: p 2 }) exp (-S^i] - S 1 ^]) (76) 

J geometry 1 J geometry 2 

where represents the action for the CFT (with fields ip) and the delta-functional in the 
measure represents a gluing condition between a geometry ‘1’ and a geometry ‘2’ along a 
‘cut’ which is the location, at a particular time r, of the spatial interval A : a G (—1/2,1/2) 
FI For Z ss , both geometries are that of a strip of the Euclidean plane described by complex 
coordinates ( w,w ) = a ± ir defined by boundaries at r = ±/3/4 with boundary conditions 
determined by the boundary state | Bd) introduced in (J^J) . For Z cc , both geometries are 
that of a cylinder cut of the Euclidean plane with identified boundaries at r = —/3/4, 3/3/4. 
The geometries for both Z ss and Z cc are familiar from calculations of Entanglement Renyi 
entropy (of order 2) and can be calculated from appropriate correlation functions of twist 
fields [32] which exchange two identical geometries. For Z SC1 the two glued geometries are 
different (that of a strip and a cylinder), hence the method of twist operators do not apply 
in a straightforward fashion. (See Figure [4]). In this paper, we will therefore, employ the 
method of the short interval expansion. 



Figure 4: Two different geometries, the strip and the cylinder, glued along the cut as described in the text. 
The method of the short interval expansion allows us to compute the functional integral over this geometry 
by replacing a small tube enclosing the two glued cuts by a complete basis of operators 4>ki ® (f>k 2 where the 
operators live in the two Hilbert spaces. 

The idea of the short interval expansion [33) is as follows. To begin, we express the 

20 To be precise, S[F] = 5(ipi (T<) — y> 2 (d>)) 5(^i(j4>) — ip 2 (d<)), where T< (A>) represents the limiting 
value from below (above) the cut. 
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functional integral (1751) as an overlap of two wavefunctions in IT <g> H 2 , as follows 


Z \2 = Tr(p Ai ip Aj2 ) = {lpout\lpin) = iDfp^Wi) D<p 2 (w 2 ) V£ut[<Pl> ¥> 2 ] 


' w\€Vi 


' W2 GT>2 


= Dipi(wi) /Dy9 2 M<%i|dPi - ^Pi)^Mdv 2 - <^2)(5(F[^i,^ 2 ])exp (-SVi] - S^D 
J Wl&T>i J W2&T> 2 

Vw[ 7T,^2] = /D(pi(u;i)/D^ 2 (M; 2 )5(v9i| a p 1 - Tp^)6{(f 2 \av 2 ~ ¥ 2 ) exp (-S^i] - S 1 ^]) (77) 


' wi^T>i 


' W2^T>2 


Here T> 1 (respectively, T> 2 ) is a small disc drawn around the cut in geometry 1 (respectively, 
geometry 2). 

Note that only |Tn) depends on the gluing condition since the delta functional in the 
measure does not affect |Vw)- The basic point of the short interval is that in the limit 
when the length l of the cut is small compared with the characterizing length scale of the 
geometries (in our case, when l -C (3), the wavefunction i/)i n [<pi, </? 2 ] becomes jointly localized 
at the centre (wi,wi) of the disc T>i and at the centre (w 2 ,w 2 ) of the disc V 2 Ej], and hence 
can be expanded in terms of local operators, as follows 


| An) = V, C klM <pk 1 (wi,W 1 )<pk 2 (w2,W2)\0)i®\0)2 

kl,k 2 


(78) 


Here k 2 label a complete basis of quasiprimary operators of the CFT Hilbert space. Each 
term in the sum represents a factorized wavefunction (between geometries 1 and 2), which, 
therefore, gives 1^1 


Zsc ^ ^ C'ki ,k 2 (0fci (t^l j ^l) ) str (0fc 2 (^2 j ^ 2 ) )cyl 1 
ki,k 2 

Zss ^ ^ C’ki ,k 2 (0fci (t^l, t^l) ) str (0fc 2 (^2j ^ 2 ) ) str j 

ki,k 2 

Zee ^ ^ C'ki ,k 2 (0fci (t^l, t^l) )cyl (0fc 2 (^2 j ^ 2 ) )cyl (79) 

k\ ,k 2 


Here the subscripts str and cyl refer to “strip”, and “cylinder” respectively. The one-point 
functions are evaluated on the respective geometries without any cut (see Section [2] for more 
details). The glued functional integral (1751) . (1771) is recovered by summing over ki,k 2 with 
the coefficients C7, ,k 2 ; , as clear from (1751) these are determined by the gluing condition and 
depend on the size of the cut [33J (see Section [Rll for more details). 


B.l The coefficients Ck lt k 2 

As explained in [33] (see also Section 1B|) . the coefficients Ck lt k 2 are determined by the equation 

Cfci,k 2 = hm (ziz 2 ) 2{h ^ +h ^\z 1 z 2 ) 2 ^ + ^\(j) kl (zi,zi)(j) k2 (z 2 ,z 2 ))c 2 (80) 

21—>-001,22—>-002 

21 We will take the centre of the disc in each geometry to coincide with the centre of the cut, which has 
coordinates w = ir,w = — it. 

22 In case geometries 1 and 2 are identical, the superscripts in Wi,uii,i = 1,2 indicate which sheet we are 
considering. 
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where C 2 represents two infinite planes glued along a cut A, Z 2 is the functional integral 
such a glued geometry and Z 1 is the functional integral over a single plane. This equation 
can be easily proved by inserting quasiprimary a operator at infinity in each plane in an 
equation like d76|i or (ITT)) . The two point function in the glued geometry is to be determined 
by using the uniformizing map: 


v = vTT/TWTTT 


(81) 


The normalization constants n k are determined by the following orthogonality condition of 
the quasiprimary operators 


(<l)k 1 (Zl,Z 1 )(j) k2 (z2,Z2))c = 


^ki^ki,k2 
Z 12 Z 12 


(82) 


where n kl is a normalization constant. Note that Ck lt k 2 — Ck 2 ,k 1 - Below we will use the 
notation 


Ck 1} k 2 — Cfc 1 ,fc 2 /C , o,o (83) 

Case ffci, k 2 ) = (0, 0) : We will denote the identity operator as </>o = 1. It is obvious that 

C 0 ,o = Z 2 /Zl (84) 


Case (fci ,k 2 ) = (k, 0 ) : The only case where Ck ,0 7 ^ 0 is when (faXz, z ) is a descendent of the 
identity operator, e.g. T(z), T{z ), : T(z)T(z) :, A(z), A(z) etco E.g. 


Ct,o — C T ,o/C 0}0 — C T> 0 — —; C TT0 — .. 


(85) 


All other Ck ,0 vanish as they are proportional to a one-point, function of a primary operator 
on the Riemann surface (and hence to that on the complex plane). 

Case (Zi.Z?) = (primary, primary) : In case 4 > k 1 ,4>k 2 are primary operators, (l80|) gives 


1 /Ze* 7 r/2 \ 

Cfci k 2 Cfci,fc 2 


( 86 ) 


Case (Zi.fc?) = (descendent, descendent) : Incase^ is of the form Z_ ni Z_„ 2 ...Z_ mi Z_ m 2 ...0z 1 
and 4>k 2 is of the form L_ ri L_ r . 2 ...Z_< il Z_< j2 ...Z; 2 , we can show that 


C kl ,ha = m ,£ s A(ni, n 2 , ■~,rn 1 ,m 2 , ...;r 1 ,r 2 , ...,si,s 2 , ...)l 2(hk i +hk i\ 

h kl = h h + n, h k2 = h h + ^ m (87) 


where A(...) is a numerical coefficient. 

23 Here A (z) = :TT \{z) — ^ d\T is the level 4 quasiprimary descendent of the identity. 
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